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Abstract 

in 

In many physical applications solitons propagate on supports whose topo- 
logical properties may induce new and interesting effects. In this paper, we 
investigate the propagation of solitons on chains with a topological inhomo- 
\ geneity generated by the insertion of a finite discrete network on the chain. 

o 

' For networks connected by a link to a single site of the chain, we derive a gen- 

JO ' eral criterion yielding the momenta for perfect reflection and transmission of 

traveling solitons and we discuss solitonic motion on chains with topological 



inhomogeneities . 

> 

In the last decades, a huge amount of work has been devoted to the study of 
the propagation of discrete solitons in regular, translational invariant lattices. 
However, in several systems, like networks of nonlinear waveguide arrays, Bose- 
Einstein condensates in optical lattices, arrays of superconducting Josephson 
junctions and silicon-based photonic crystals, one can engineer the shape (i.e., 
the topology) of the network. Correspondingly, an interesting task is the study 
of the propagation of solitons in inhomogeneous networks. The general idea 
of this work is that network topology strongly affects the soliton propagation. 
We provide a general argument giving the momenta of perfect transmission and 
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reflection for a soliton scattering through a finite general network attached to a 
site of a chain: the momenta of perfect transmission and reflection are related 
in a simple way to the energy levels of the attached network. This criterion 
directly links the transmission coefficients with the network adjacency matrix, 
which encodes all the relevant informations on its topology. Such relation puts 
into evidence the topological effects on the soliton propagation. The situations 
where finite linear chains, Cayley trees and other simple structures are attached 
to a site of an unbranched chain are investigated in detail. 

I. INTRODUCTION 

The analysis of nonlinear models on regular lattices [1-3], as well as the investigation 
of linear models on inhomogeneous and fractal networks [4] has attracted a great deal of 
attention in the last decades: while nonlinearity dramatically modifies the dynamics, allow- 
ing for soliton propagation, energy localization, and the existence of discrete breathers [2], 
topology mainly affects the energy spectrum giving rise to interesting phenomena such as 
anomalous diffusion, localized states, and fracton dynamics [4]. It is now both timely and 
highly desirable to begin a thorough investigation of nonlinear models on general inhomo- 
geneous networks, since one expects not only interesting new phenomena arising from the 
interplay between nonlinearity and topology, but also an high potential impact for appli- 
cations to biology [5,6] and to signal propagation in optical waveguides [7]. Recently, the 
effects of uniformity break on soliton propagation [8,9] and localized modes [10] has been 
investigated by considering Y-j unctions [8,9] (consisting of a long chain inserted on a site of 
a chain yielding a star-like geometry) or geometries like junctions of two infinite waveguides 
or the waveguide coupler [10]. Here, we consider general finite networks inserted on a chain. 

The discrete nonlinear Schrodinger equation (DNLSE) is a paradigmatic example of 
a nonlinear equation on a lattice which has been successfully applied to several contexts 
[11,12]: in particular, it has been used to describe the physics of arrays of coupled optical 
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waveguides [13,14] and arrays of Bose-Einstein condensates [15]. It is well known that, on a 
homogeneous chain, the DNLSE is not integrable [12]; nevertheless, soliton-like wavepackets 
can propagate for a long time and the stability conditions of soliton-like solutions can be 
derived within variational approaches [16]. Furthermore, the dynamics of traveling pulses 
has been investigated in detail in literature [17-19] (more references are in [11]). The simplest 
example of an inhomogeneous chain is provided by an external potential localized on a site of 
the chain [20-23,11]: an experimental set up with a single defect has been recently realized 
with coupled optical waveguides [24]. Another relevant example of an inhomogeneous chain 
is obtained adding an additional Fano degree of freedom coupled to the a site of the chain, 
which gives the so-called Fano-Anderson model (see [25-27] and references therein). 

As a first step in the investigation of the properties of nonlinear models on general 
inhomogeneous networks, we shall analyze the propagation of DNLSE solitons on a class of 
inhomogeneous networks built by suitably adding a finite topological inhomogeneity to an 
unbranched chain. The general framework where this analysis can be carried out is provided 
by graph theory [28]; in particular, we shall consider networks where a finite discrete graph 
G° is attached by a link to a site of the homogeneous, unbranched chain (see Fig.l) while 
all the sites potentials q are set to a constant. Such systems may be experimentally realized 
by placing the nonlinear waveguides in a suitable inhomogeneous arrangement, like the one 
depicted in Fig.l. We mention that in arrays of Bose-Einstein condensates one can build 
up geometries, which differ from the unbranched chain, by properly superimposing the laser 
beams creating the optical lattices [29]. In superconducting Josephson arrays [30], present- 
day technologies allow for to prepare the insulating support for the junctions in order to 
create structures of the form "chain + a topological defect". In the context of coupled 
nonlinear waveguides [3], one should couple the waveguides according the geometry of the 
graph G°, and couple this network of waveguides to a single waveguide of the unbranched 
chain; a similar engineering should be requested to realize photonic crystal circuits [31] 
obtained merging the circuit G° to the unbranched chain. 

As we shall discuss in Section IV, the shape of the attached graph G° affects the trans- 
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mission and reflection coefficients as a function of the soliton momentum: as an example 
of this general phenomenon, we consider unbranched chains to which simple graphs, like 
finite chains and Cayley trees, are added. Our analysis points to the fact that the topology 
of the network (i.e., how its sites are connected) controls the transmission properties, and 
that one can modify soliton propagation by varying the topology of the inserted network. 
In particular, we shall show that the momenta of perfect transmission are determined by 
the energy levels of the inserted graph, i.e., by the eigenfrequencies of the G°'s oscillation 
modes in the linear case (see Section IV). 

On a chain, stable solitonic wavepackets can propagate for long times [11]. When a graph 
G° is inserted, one can study how the presence of this topological inhomogeneity modifies 
the soliton propagation. We numerically evaluate the transmission coefficients and we com- 
pare the numerical results with analytical findings obtained for a relevant soliton class, to 
which we refer as large-fast solitons. For large fast solitons the transmission coefficients 
can be evaluated within a linear approximation [25]. Indeed the characteristic time for the 
soliton-topological defect collision are very small with respect to the soliton dispersion time; 
therefore, the soliton scattering can be approximated by the scattering of a plane-wave. 
However, as we numerically checked in the Figs. 2-6, the nonlinearity still plays a role, giving 
long-lived solitons, especially near the momenta of perfect reflection or perfect transmission: 
it keeps the soliton shape during its propagation. Since in many experimental settings one 
can easily check if the reflected wavepacket is vanishing, this work could provide a basis for 
a topological engineering of solitonic propagation on inhomogeneous networks. 

The plan of the paper is the following. In Section II we review the properties of the 
DNLSE on a chain and we introduce the variational approach to investigate the soliton 
dynamics. In Section III, by using graph theory [28], we define the DNLSE on inhomogeneous 
networks built by adding a topological perturbation to an unbranched chain; furthermore, we 
explain the numerical techniques used in the paper for the study of the soliton scattering, 
and we discuss the range of validity of the linear approximation used in the analytical 
computations. In Section IV we present our analysis yielding the conditions on the spectrum 
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of the finite graph G° in order to obtain total reflection and transmission of solitons. In 
Section V we study the relevant case when G° is a finite, linear chain and we show that the 
Fano-Anderson model [25,26] can be realized within our approach by considering a single 
link attached to the unbranched chain. In Section VI we show that, for self-similar graphs 
G°, the values of momenta for which perfect reflection occurs becomes perfect transmission 
momenta when the next generation of the graph is considered and we study in detail the 
case of Cayley trees as an example of self-similar structures. In Section VII we study the 
transmission coefficients for three different inserted finite graphs: loops, stars and complete 
graphs. Finally, Section VIII is devoted to our concluding remarks. 



II. DNLSE ON A CHAIN 

Besides its theoretical interest, the DNLSE describes the properties of interesting sys- 
tems, such as arrays of coupled optical waveguides and arrays of Bose-Einstein condensates. 
On a chain the DNLSE reads 

dijj n 1 

i-j^r = --(il>„+i + i> n -i) + A | i>„ | 2 + e n ij n (l) 

where n is an integer index denoting the site position and the normalization condition is 
J2 n I VVi | 2= 1- 111 Eq.(l) one has a kinetic coupling term only between nearest-neighbour 
sites, but the effect of next- nearest-neighbour coupling and long-term coupling has been 
also often considered (see the reviews [3,11] for more references): in the present paper we 
consider only a constant nearest-neighbour interaction, but from the next Section we allow 
for that the number of nearest-neighbours of a site is not constant across the network (like 
for the simple chain), but it can vary according the topology of the graph. 

In condensate arrays, ip n (j) is the wavefunction of the condensate in the nth well. Time 
r is in units of Ti/ZK, where K is the tunneling rate between neighbouring condensates; 
e„ = E n /2K where E n is an external on-site field superimposed to the optical lattice and 
the nonlinear coefficient is A = U /2K, where U is due to the interatomic interaction and it 
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is proportional to the scattering length (U is positive for 87 Rb atoms and is negative for 7 Li 
atoms) . 

In arrays of one-dimensional coupled optical waveguides [13] ■0n( r ) is the electric field in 
the nth-waveguide at the position r and the DNLSE describes the spatial evolution of the 
field. The parameter A is proportional to the Kerr nonlinearity and the on-site potentials 
e n are the effective refraction indices of the individual waveguides. As the light propagates 
along the array, the coupling induces an exchange of power among the single waveguides. In 
the low power limit (i.e. when the nonlinearity is negligible), the optical field spreads over 
the whole array. Upon increasing the power, the output field narrows until it is localized 
in a few waveguides, and discrete solitons can finally be observed [13,14]. Experiments 
with defects (i.e., with particular waveguides different from the others) have been already 
reported [24]. 

On a chain DNLSE soliton-like wavepackets can propagate for a long time even if the 
equation is not integrable [3]. Let us consider, at r = 0, a gaussian wavepacket centered in 
£(t = 0) = £oi with initial momentum k and width 7(7- — 0) = 7 : its time dynamics are 
studied resorting to the Dirac time-dependent variational approach [32] which well repro- 
duces the exact results in the continuum theory [33] . In its discrete version the wavefunction 
can be written as a generalized gaussian 

Mr) = Vie • e -^**(«-«*§(»-« 2 (2) 

where £(r) and 7(7-) are, respectively, the center and the width of the density p n =\ ip n | 2 , 
and k{r) and 5{r) are the momenta conjugate to £(r) and 7(7-) respectively; K, is just a 
normalization factor. The wave packet dynamical evolution is obtained from the Lagrangian 
C = J2 n i^n^n ~ with the equations of motion for the variational parameters £, 7, k, 5. 
In the absence of external potential (e n = 0), one obtains the Lagrangian [34] 



A/C 



C = K, £ e -( 2 " + "« +2 «- 2 « +1 ^ cos[5(n+l/2-0 + A;]-^- £ e 



4(n-£) 2 /7 2 

2 



n=— 00 



00 ( i \ 

+ IC E ~2 {n " ° 2 + " " k(n " + *M e_2(n_?)2/72 - ( 3 ) 



With 7 not too small (7 ^> 1), we can replace the sums over n with integrals: to evaluate 
the error committed, we recall that [35] 



E e 

o = "°° _ (B _ oa =1 + Q(e-^). (4) 



J dn e 



7 



In this limit the normalization factor becomes fC = y2/7T7 2 . We finally get [15] 

where n = I/27 2 + 7 2 5 2 /8. The equations of motion are 

k = (6) 

£ = sink-e- v (7) 

5 = cos fc (4/ 7 4 - <5 2 ) e~ f ' + 2A/ v^F 7 3 (8) 

7 = 7<5cos fc • e~ v : (9) 

fc(r) = k is conserved. Notice that, due to the discreteness, the group velocity cannot 
be arbitrarily large (£ sin A; < 1). As Eq.(4) clearly shows, the variational equations 
of motions (9) are meaningful only for large solitons: the Peierls-Nabarro potential does 
not appear in (9), and the equations feature momentum conservation. We mention that in 
uniform DNLSE chains a threshold condition for the soliton propagation appears: only if the 
soliton is sufficiently broad solitons may freely move [36]. In the following we shall consider 
only large-fast solitons, so that the variational equations of motions (9) are appropriate; 
however, to study the propagation of localized discrete breathers in inhomogeneous networks 
one should study the Lagrangian (3). 

When 7 = and 5 = 0, the shape of the wavefunction does not vary and one has 
a variational soliton-like solution where the center of mass move with a constant velocity 
£ = constant. If A > the conditions 7 = and 5 = can be satisfied only if cos k < 
(i.e., only when the effective mass is negative): for this reason in the following we take only 
momenta n/2 < k < n (positive velocities) or — n < k < —n/2 (negative velocities). In 



particular, for 5(r = 0) = 5 = and large enough solitons (70 3> 1), the condition on A 
allowing for a soliton solution is [15] 

~ 1 — cos k . _ x 

A soi « 2^ L (10) 

7o 

The stability of variational solutions has been numerically checked showing that the 
shape of the solitons is preserved for long times. In the following, we use the term "solitons" 
to name the solutions of the variational equations (6)- (9). One can have a similar criterion 
using other variational approaches [16]. One also expects that the integrable version of the 
DNLSE, the so-called Ablowitz-Ladik equation [37,12], provides results very similar to those 
obtained in this paper. 



III. DNLSE ON GRAPHS 

The DNLSE (1) can be generalized to a general discrete network by means of graph 
theory. A graph G is given by a set of sites % connected pairwise by set of unoriented links 
(i, j) defining a neighbouring relation between the sites. The topology of a graph is described 
by its adjacency matrix A it j which is defined to be 1 if i and j are nearest-neighbours, and 
otherwise. The DNLSE on a graph reads as 

i-Q^ = - g S Ai ^i + A I ^ | 2 ^ + erfi (11) 

Equation (11) describes the wavefunction dynamics in a wide range of discrete physical 
systems, and it can be applied to regular lattices as well as to inhomogeneous networks 
such as fractals, complex biological structures and glasses. The properties of Eq.(ll) on 
small graphs has been investigated in [38]. The first term in the right-hand side of Eq.(ll) 
represents the hopping between nearest-neighbours (with tunneling rate proportional to 
Aij), the second is the nonlinear term, and the third one describes superimposed external 
potentials. We remark that in Eq.(ll) the numbers of nearest neighbours is site-dependent. 
Furthermore, since A it j — 1 if i and j are nearest-neighbours, one is assuming that the 
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tunneling rate between neighbouring sites entering Eq.(ll) is constant across the array and 
it is (in the chosen units) equal to 1. Below, we shall consider also the case of a tunneling 
rate which is not constant across the array. 

We shall focus on the situation where the graph G is obtained by attaching a finite graph 
G° to a single site of the unbranched chain (see Fig.l) and setting = for all the sites. 
We denote the sites of the unbranched chain and of the graph G° with latin indices m,n, . . . 
and greek indices a, (3, • • • respectively. A single link connects the site n = of the chain 
with the site a of the graph G°. 

The scattering of a soliton through the topological perturbation can be numerically stud- 
ied as follows. At t = (hereafter, we refer to r as a time even if for the optical waveguides 
it represents a spatial variable) one prepares a gaussian soliton (2) centered well to the left 
of (i.e., £ < 0) moving towards n = (sin(/c) > 0) with a width related to the nonlin- 
ear coefficient according to Eq.(10). From Eq.(ll) the time evolution of the wavefunction 
may be numerically evaluated: when t s £ / sin(/c) the soliton scatters through the finite 
graph G° (sm(k) being the group velocity of the soliton). At a time r well after the soliton 
scattering (i.e. r t s ), the reflection and transmission coefficients 1Z and T are given by 

ft = £ I Mr) I 2 (12) 

n<0 

T=Y,\Mr)\ 2 - (13) 

«>o 

Note that, while in the linear case (A = 0) one has 1Z + T — 1, in general, nonlinearity 
violates unitarity by allowing for phenomena such as soliton trapping; nevertheless, there 
are regimes where soliton trapping is negligible and 1Z + T 1. We numerically checked 
that in the time dynamics reported in the paper this condition is well satisfied. Situations 
corresponding to resonant scattering (i.e. TZ = or T = 0) have a particular relevance: 
in fact, these situations can be easily experimentally detected, and the soliton-like solution 
is stable also well after the scattering, as it is numerically verified in different examples of 
resonant reflection and transmission. 
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For an important class of soliton solutions (to which we refer as large-fast solitons) the 
scattering through a topological inhomogeneity can be analytically studied using a linear 
approximation. The interaction between the soliton and the topological inhomogeneity is 
characterized by two time-scales: the time of the soliton-defect interaction r int = 7/ sin A; 
and the soliton dispersion time (i.e. the time scale in which the wavepacket will spread 
in absence of interaction) Tdi sp = 7/(4sin (I/27) cosA;) [25]. For large (7 ^> 1, as in many 
relevant experimental settings) and fast solitons, i.e. 

v = sin k > (2/7) cos k, (14) 

one has that the soliton may be considered as a set of non interacting plane waves while 
experiencing scattering on the graph; thus the soliton transmission may be studied by con- 
sidering, in the linear regime (i.e., A = 0), the transport coefficients of a plane wave across 
the topological defect. The use of the linear approximation for the analysis of the interac- 
tion of a fast soliton with a local defect in the continuous nonlinear Schrodinger equation 
is reported in [39]. Later, we shall compare the analytical findings with a numerical solu- 
tion of Eq.(ll), namely with the reflection and transmission coefficients 1Z and T given by 
Eqs.(12)-(13). 

IV. A GENERAL ARGUMENT FOR RESONANT TRANSMISSION 

In this Section we show that, if a large- fast soliton scatters through a topological pertur- 
bation of an unbranched chain, the soliton momenta for perfect reflection and transmission 
are completely determined by the spectral properties of the attached graph G°: in particular 
one has 1Z — 1 if 2 • cos A; coincides with an energy level of G°, while T = 1 if 2 • cos A; is 
an energy level of the reduced graph G r , i.e., of the graph obtained from G° by cutting the 
site a from G° (see Fig.l). In algebraic graph theory (see e.g. [28,40]) the energy level of a 
graph is simply defined as an eigenvalue of its adjacency matrix. We will call A and A r the 
adjacency matrices of G° and G r , respectively. 
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For large-fast solitons, the pertinent eigenvalue equation to investigate is 

-oEAA'^rt (15) 

(here % is a generic site of the network) . The solution corresponding to a plane wave coming 
from the left of the chain is ip n = ae' lkn + be~ lkn for n < and ip n = ce lkn for n > 0, so 
that pi = — cos A;. The reflection coefficient is given by 7Z —\ b/a | 2 and the transmission 
coefficient by T =| c/a | 2 . The continuity at requires a + b = c. The equation in is 

-hae~ ik + be lk + ce ik + ip a ) = - cos k ■ (a + b) (16) 

while in a one has 

~ha + b+ E A °«M = -cosfc-^ Q . (17) 

z v eG0 

At the sites 77 of G r one obtains 

E 4UVy = -cosfc.^. (is) 

z v 'eG° 

For perfect reflection, i.e. for the momenta fc's such that 7£(/c) = 0, one has c = 0, 
a = —b and from Eq.(16) = —2a sin k. Therefore Eqs.(17) and (18) reduce to the 
eigenvalue equation for the adjacency matrix A and, apart from the trivial case cos(A;) = 0, 
they are satisfied only if 2 cos A; coincides with an eigenvalue of A . 

At variance, in order to find the momenta fc's such that T(k) = (perfect transmission), 
one has b = 0, a = c and from Eq.(16) i\) a = 0. Therefore, Eqs.(18) reduces to the eigenvalue 
equation for A r and it is satisfied only if 2 cos A; coincides with an eigenvalue of A r . 

This general argument can be easily extended to the situation where p identical graphs 
G° are attached to n = 0: indeed, now one has only to replace in Eq. (16) ip a with pip a and 
the conditions for T(k) = and 71(h) = do not change. 

The stated result holds for the case where the tunneling rates between the neighbour sites 
of G° are constant (and equal to 1 in the chosen units): however one can also consider in G° 
non-uniform tunneling rates t® , > 0, where 77 and rj' are nearest-neighbour sites belonging 
to G° (t° , = if Al. = 0). The DNLSE at a site ri of G° becomes 
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^ = -2£<WV+ A l^| 2 ^ (19) 
if 

while the DNLSE at the sites n of the of the chain remain unchanged. Now, the criterion 
states that 1Z = 1 if 2 cos A; coincides with an eigenvalue of the matrix i° ,, while T = 1 if 
2 cos is an eigenvalue of the matrix ii „, defined as tl „, = t°„ if r> and 7/ belong to the 
reduced graph G r . 



V. FINITE LINEAR CHAINS 

As a first simple application of the argument given in Section IV, we consider a single 
site a attached via a single link to the site of the unbranched chain. For Bose-Einstein 
condensates in optical lattices [29] the setup "infinite chain + single link" or "infinite chain 
+ finite chain" may be realized by using two pairs of counterpropagating laser beams to 
create a star-shaped geometry in the x-y plan [41] and manipulating the frequencies of the 
superimposed harmonic magnetic potential so that in the y direction only few sites can be 
occupied. In an optic context, chains of coupled waveguides are routinely built and studied 
[7]: one can obtain the configuration "infinite chain + single link" by coupling a further 
waveguide to a waveguide of the chain. 

The DNLSE in the sites and a reads 

i^T = ~\(^ + ^-1 + ^a) + A I ^0 I 2 ^0 (20) 

and 

i^ = ~lfo + A|iM 2 ifc,. (21) 

It is transparent from Eqs.(20)-(21) that the wavefunction ip a (r) may be interpreted as 
an additional local Fano degree of freedom, yielding the so-called Fano-Anderson model 
[25-27]. In our approach, such degree of freedom is interpreted as a single link attached to 
the unbranched chain. As it is well known, the Fano-Anderson model describes interesting 
scattering properties: adding a generic finite graph (instead of a single link) gives rise to 
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a yet richer variety of behaviors. We mention that in [27] the Fano degree of freedom is 
coupled to several sites of the unbranched chain: this would correspond in our description 
to a site linked to several sites of the chain, and, in general, to graphs attached to several 
sites of the chain. For simplicity, in the following we limit ourself to graphs inserted in a 
single site of the unbranched chain. 



For large-fast solitons, when a single link is added to the unbranched chain, the reflection 
coefficient K from Eqs.(20)-(21) is found to be [25] 



in the regime where Tdi sp ^> T int , we verified that the numerical results of the soliton scat- 
tering against the link are in agreement with Eq.(22) (see Fig. 2). One sees that, when k is 
approaching n (solitons becoming slower), the agreement becomes worse. From the general 
results, there are no fully transmitted momenta and 1Z — 1 only for k = n/2 and k = n, as 
one can also see by a direct inspection of (22). 

One can also attach a finite chain of length L at the site 0. In the linear approximation 
for large-fast solitons (i.e., Eq.(15)), the solution corresponding to a plane wave coming 
from the left of the unbranched chain is tp n = ae lkn + be~' lkn for n < and ipn = ce lkn for 
n > 0, while in the attached chain ip a = fe lka + ge~ tka (a = 1, • • • , L denotes the sites of 
the attached chain, and a — 1 is the site linked to n — 0). Eq.(15) for n = 0, a = 1, and 
a = L — 1 yields respectively 



1 



(22) 



K = 



l + 4sin 2 (2A;) ' 



a + b = c = f + g 



(23) 



-{ae~ lk + be~ ik + ce lk + fe lk + g~ lk ) = fi{a + b) 



(24) 



and 



h fe ik(L-l) + ge -ik(L-l)j = ^ fe ikL + ge -ikL^ 



(25) 



where n = — cos k. 
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We have five unknowns (a, b, c, f, and g) and four equations (23)- (24) (the remaining 
condition being provided by the normalization). One may easily determine b/a, c/a, f /a, 
and g/a, getting 

b e Hk( e 2ikL _ ^ 



a 1 — 2e 2ik + e 2 * fc ( L + 2 ) ' 

which leads to 



(26) 



= l /a 1 = [cos (kL) - cos {k{L + 2))] 2 + sin 2 (kL) ' ( ^ 

and T =| c/a | 2 = 1-7?.. For L — 1, Eq.(27) reduces to Eq.(22). In agreement with 
the general argument of Section IV, the number of minima and maxima increases with L. 
Eq.(27) for L = 2 is compared in Fig.2 with the numerical results. 

We notice that in the limit L — > oo the considered problem corresponds to the propa- 
gation of a soliton in a the so-called star graph, which has been recently investigated in the 
context of two-dimensional networks of nonlinear waveguide arrays [8] and Y-j unctions for 
matter waves [9]. 



VI. CAYLEY TREES 

Eq.(27) yields that the values of k allowing for perfect reflection (TZ = 1) for the length 
L coincide with the momenta of full transmission (T = 1) when G° is a chain of length 
L + 1. This property is readily understood: in fact, if G° is a chain of length L the fc's for 
which 1Z = 1 correspond to the energy levels of G°. If G° is a chain of length L + 1, the 
values of fc's for which T = 1 correspond to energy levels of the reduced graph G r , which, 
in this case, is again given by a chain of length L. This is clearly a general property of any 
self-similar graph. As an example, we study in this section the situation in which G° is a 
Cayley tree of branching rate p and generation L (see Fig. 3). 

Let us consider the linear approximation for large-fast solitons Eq.(15). The plane wave 
coming from the left of the unbranched chain is ip n = ae lkn + be~ lkn for n < and ip n = ce lkn 
for n > 0. This fixes /i = — cos k. Furthermore, from the continuity in it follows a + b = c. 
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For a Cayley tree, the eigenfunction must have, by symmetry, the same value at all the 
sites belonging to the same generation. If we denote by ipp the eigenfunction at the sites at 
distance j3 — 1, • • • , L from n = 0, the eigenvalue equation (15) at the site a — 2, • • • , L — 1 
reads 

~(^a-l + WWl) = V>lf>a- (28) 

The plane wave solutions of Eq.(28) can be written as 

^ = ^(fe tk ' a + ge- Wa ), (29) 

where a = 1, • • • , L: in this way one gets from Eq.(28) /i = —y/pcosk', so that k' = 
arccos (p _1//2 cos A;). Eq.(15) for n = 0, a = 1, and a = L — 1 gives respectively 

-Uae~ lk + be ik + ce ik ) - ^-(fe* k ' + ge~ ik ') = fi(a + b) (30) 
2 

~\{a + b)- l -{fe 2lk ' + ge- 2M ) = ^{fe M + ge^) (31) 

and 

_I (/e i*'(L-i) + ge -ik'(L-i)^ = JL( fe ik>L + ^-ifc'i). (32) 

2 y/P 

Using Eqs.(30)-(32) and the condition a + b = c, one can determine b/a, c/a, f /a, and g/a: 
the resulting expressions is rather involved and here we will not explicitly write them. In 
Fig. 3 we plot the numerical and analytical results for the reflection coefficient TZ when two 
Cayley trees, respectively with L = 5 and L = 6, are attached to the unbranched chain. 
One sees that when one pass to the next generation, the momenta for which full reflection 
occurs becomes momenta of full transmission. 



VII. FURTHER EXAMPLES OF INSERTED GRAPHS 



In general, it is possible to consider the scattering of a large-fast soliton through a large 
variety of inhomogeneous networks. Here, we analyze three further examples of network 
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topologies: stars, loops and complete graphs. In each case, the reflection coefficients for 
large-fast solitons are derived with a procedure analogous to the one adopted in the previous 
section for Cayley trees. The analytical findings are compared with numerical results. 

Loops: Let us consider a loop graph, i.e., a finite chain of L sites oti, • • • , oll such that 
the sites «i and are linked to the site n = of the unbranched chain. For large-fast 
solitons, the reflection coefficient TZ in the linear approximated regime is 

n = 2 [i + C Q S (A;(L-l))] 2 + sin 2 (A;(L-l)) 

6 — 2 cos (2k) + 5 cos k cos (kL) — cos (k(L + 3)) + sin k sin (kL) ' 

We note that the number of momenta of perfect reflection and perfect transmission increases 
with L. Furthermore, at large L, the transmission properties of the loops become similar 
to those of a finite (not closed) chain [see Eq.(27)]. In Fig. 4 the numerical and analytical 
results are compared and a figure of the loop graph with L = 3 is provided. 

Stars: The p-star is the graph composed by a central site linked to p sites, which are in 
turn connected only to the central site. Let us consider the case where G° is a p-star graph 
and a is the center, linked to n — 0. In the linear approximation we have 

K = ? . (34) 

l+(^-(P-Dtan^) 2 

In Fig. 5 the numerical and analytical results are compared. Eq.(34) shows that perfect 
transmission (TZ = 0) is obtained only for for the momentum k = tt/2. This can be directly 
proved applying the criterion of Section IV: indeed G° is a p-star, while G r consists of p 
disconnected sites. Therefore, all the eigenvalues of the adjacency matrix A r equal zero, and 
the only momentum of perfect transmission is k = ir/2. 

Complete graphs: The complete graph Km of M sites is the graph where every pair of 
sites in linked [28]: e.g, K 3 is a triangle. Inserting K L+ i at the site n = of the unbranched 
chain (so that n = is one of the sites of K L+1 ), one gets: 

n = i+4 (L - i - L 2 r fc)2 shT? (35) 
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For L » 1, 7£ « 1/(1 + 4 sin 2 fc), therefore the complete graph behaves as a single effective 
defect [compare with Eq.(22)]. The comparison between the numerical and analytical results 
is presented in Fig. 6. 

VIII. CONCLUSIONS 

As a first step in addressing the issue of the interplay between nonlinearity and topology 
we studied the discrete nonlinear Schrodinger equation on a network built by attaching 
to a site of an unbranched chain a topological perturbation G°. The relevant situation 
corresponding to the Fano-Anderson model is obtained when one considers a single link 
attached to the linear chain. We showed that, by properly selecting the attached graph, 
one is able to control the perfect reflection and transmission of traveling solitons. We 
derived a general criterion yielding - once the energy levels of the graph G° is known - the 
momenta at which the soliton is fully reflected or fully transmitted. For self-similar graphs 
G°, we found that the values of momenta for which perfect reflection occurs become perfect 
transmission momenta when the next generation of the graph is considered. For finite linear 
chains, loops, stars and complete graphs, we studied the transmission coefficients and we 
compared numerical results form the discrete nonlinear Schrodinger equation with analytical 
estimates. Our results evidence the remarkable influence of topology on nonlinear dynamics 
and are amenable to interesting applications in optics since one may think of engineering 
inhomogeneous chains acting as a filter for the motion of soliton [42] . 

Acknowledgments: We thank M. J. Ablowitz, P. G. Kevrekidis and B. A. Malomed for 
discussions. 
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FIGURES 




FIG. 1. Inserting a graph G on a site of a linear chain: the points of the chain are denoted with 
integers n and the point in which the graph is attached is n = 0; a is the point of G° connected to 
0. G r is obtained subtracting a from G°. 
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FIG. 2. Reflection coefficient 1Z as a function of k (with k between tt/2 and it) when a chain 
with length 1 (i.e., a single link) and 2 are attached. Empty circles (L = 1) and stars (L = 2) 
correspond to the numerical solution of Eq.(ll): in this figure, as well in the followings, as initial 
condition we choose a Gaussian with initial width 70 = 40 and momentum k. Solid lines correspond 
to the analytical prediction (27). 
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FIG. 3. Reflection coefficient 1Z as a function of k when a Cayley tree with length 5 and 6 
are attached. Empty circles (L = 6) and stars (L = 5) correspond to the numerical solution of 
Eq.(ll). Solid lines correspond to the analytical prediction (see text). As required from the general 
argument in Section III, the values of k for which one has perfect reflection (1Z(k) = 1) for L = 5 
correspond to perfect transmission (IZ(k) = 0) for L = 6. 




FIG. 4. Reflection coefficient 1Z as a function of k when loops with L = 4, 7 and 10 are inserted 
at a site of the unbranched chain. Empty circles (L = 4), squares (L = 7) and diamonds (L = 10) 
are obtained from the numerical solution of Eq.(ll). Solid lines correspond to the analytical 
prediction (33). The small figure represents a loop with L = 3 inserted in the unbranched chain. 
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FIG. 5. Reflection coefficient 1Z as a function of k when stars with p = 3, 4 and 6 are inserted. 
Empty circles (p = 2), squares (p = 3) and diamonds (p = 5) are obtained from the numerical 
solution of Eq.(ll). Solid lines correspond to the analytical prediction (34). The inset represents 
the situation where the attached graph G° is a star with p = 3. 
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FIG. 6. Reflection coefficient TZ as a function of k when the complete graphs K% and K\\ arc 
inserted at a site of an unbranched chain. Empty circles (K3) and squares (Kn) are obtained from 
the numerical solution of Eq.(ll). Solid lines correspond to the analytical prediction (35). 




23 



